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We present the Belavkin filtering equation for the intense balanced hetero- 
dyne detection in a unitary model of an indirect observation. The measuring 
apparatus modelled by a Bose field is initially prepared in a coherent state and 
the observed process is a diffusion one. We prove that this filtering equation 
is relaxing: any initial square-integrable function tends asymptotically to a 
coherent state with an amplitude depending on the coupling constant and the 
initial state of the apparatus. The time-development of a squeezed coherent 
state is studied and compared with the previous results obtained for the 
measuring apparatus prepared initially in the vacuum state. © 2012 Optical 
Society of America 
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1. Introduction 

The theory of quantum measurements continuous in time is firmlv based on quantum stochas- 
tic calculus (QSC) developed by Hudson and Parthasarathy EJ- The time development of 
the posterior state conditioned by a trajectory, of the results of a continuous measurement 
is given by the Belavkin filtering equation |3|46j|. The measurement is taken on a Bose field 
interacting with the quantum system in question and enables one to perform its indirect 
observation. The Bose field can be treated as an approximation to the electromagnetic field. 
The filtering equation has the form of the Ito quantum stochastic differential equation and 
plays the role analogous to that of the Schrodinger equation for an unobserved quantum 
systems. The observed process has the properties of a diffusion or/and a counting one. The 
Belavkin filtering equation was obtained under the assumption that the Bose field modelling 
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the apparatus was initially prepared in the vacuum state. In js- 10| this assumption has been 
relaxed and the filtering equation has been derived for the Bose field prepared in a coherent 
state for the counting and diffusion observations. 

The aim of this paper is twofold — we present the Belavkin equation for the diffusion 
observation with the apparatus prepared initially in a coherent state and we discuss the time 
development of a squeezed coherent state undergoing a balance heterodyne measurement 



In contrast to the result of 12], showing that for the apparatus initially prepared in 
the vacuum state a single and double heterodyne detection does not destroy the squeezed 
coherent state and drives the system asymptotically to the vacuum state, the state asymp- 
totically relaxes to the coherent one with the amplitude independent of the initial state of 
the system. Next, we generalize this observation and prove that any square-integrable ini- 
tial wave function relaxes to the coherent state with the given amplitude depending on the 
strength of the interaction between the system (single-mode field) and the apparatus and 
the initial state of the apparatus. Consequently, the considered filtering equation describes 
in fact a control of the quantum system by driving its state to the coherent one with the 
given amplitude. 

Though the methods of QSC are neither widely used nor well-known by physicists, they 
deserve in our opinion more interest as the effective tools in modelling physical systems in- 
teracting with measuring devices. We refer the readers who not familiar with the theory of 



quantum measurements continuous in time to the books [13Ml5| or the recent papers 16Ml9l|. 
More rigorous approach to the subject can be found in [20l-|23l|. The experimental achieve- 
ments in the area of continuous in time observations of quantum systems were reported, for 
instance, in 
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27|. 



The paper is organized as follows. In Section 2 we present the basic rules of quantum 
stochastic calculus. Section 3 is devoted to the presentation of a linear version of the Belavkin 
filtering equation for a balanced heterodyne detection of the diffusion type for the apparatus 
prepared initially in the coherent state. We choose the linear form of the filtering equation 
for an unnormalized posterior wave function) instead of the nonlinear version derived in 



12]], because it is more convenient to deal with. The physical interpretation of the linear 



filtering equation one can find, for instance, in [221 . |28| |. We put forward here the approach 
of the generating map for the underlined continuous observation. In Section 4 we discuss the 
time development of a coherent and a squeezed coherent state evolving under a continuous 
diffusion observation of a coherent channel. 

2. Quantum stochastic calculus 

In this section we recall some basic rules of quantum stochastic calculus (QSC) in the boson 
Fock space 0, 0| • Denote by J 7 the (symmetrical) Fock space over the Hilbert space /C = 
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C n ® L 2 (M + ) of all square integrable functions from M + into C". For any / £ JC one can 
define a coherent vector, e(f), by the formula 

e(f) = exp(~\\f\\l\ (l,/,(2!)- 1 / 2 /®/,(3!)- 1/2 /®/® /,•••) • (1) 

In particular, e(0) = (1,0,0,...) £ J 7 is the Fock vacuum. The annihilation, creation and 
number processes: Bi{t), B]{t) and A^t) are defined on the dense in J 7 linear span of all 
coherent vectors [l|, Q] as follows: 

Bj(t)e(J) = f f)(s)ds e(f) , (2) 



o 



e{g)\B]{t)e{f)) = / 9j (s)ds (e(g)\e(f)) , (3) 



o 



e{g)\A l3 (t)e(f)) = / g l (s)f J (s)ds(e(gMf)). (4) 



o 



These are the underlying processes for stochastic differential equations (QSDEs) of the Ito 
type: 

n / n \ 

dM(t) =J2(H F i*V) d Mt) + E&) dB i(t) + D,(t) dB}(t) + C{t) dt . (5) 
j=l \i=l J 

In ()5]) all the processes appearing at the Ito differentials are adapted processes on "H ® J, 
i.e. they depend on the processes (J2]|4]) up to t (present instant) and commute with the Ito 
differentials that "point to the future". If M'(t) is the process which satisfies an equation of 
the type §5$), then the differential of the product M(t)M'(t) is given by the formula 

d(M(t)M'(t)) = dM(t)M'(t) + M(t)dM'(t) + dM(t)dM'(t) . (6) 

The term dM(t) dM'(t) can be computed with the help of the multiplication table: 

d.Bi(t) dB]{t) = 6 i:j dt, dBi{t) dA kj {t) = 6 ik dBj{t), 

dA kj (t)dBj(t) = 5 ri dBl(t), dA tJ (t)dA kl (t) = 5 jk dA u (t), (7) 
and all other products vanish. 

3. Linear filtering equation for a balanced heterodyne scheme 

Let us consider a harmonic oscillator (system S) interacting with an environment modelled 
by the Bose field in a coherent state e(/). We assume that the unitary evolution operator 
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U(t ) of the compound system (system S plus one-dimensional Bose field) satisfies the QSDE 
2fJ: 



6U(t) 



y/Ji a dB ] (t) - y/Ji a) dB(t) - - a j a dt - jHdt 



U(t), U(0) 



where H = ftw (a^a + |) is the Hamiltonian of S, a is an annihilation operator, and \i G M+ 
stands for a real coupling constant. Eq. ([8]) is written in the interaction picture with respect 
to the free dynamics of the Bose field. The description of physical assumptions leading to 



this evolution can be found, for instance, in [20|, |29j]. In short, the coupling is linear in 
the field operators, the rotating-wave approximation (RWA) is made, the coupling constant 
are independent of frequency, and the spectrum of the reservoir is flat and broad. These 
assumptions are often made in quantum optics. 

Though the Bose field disturbs the free evolution of 5, it also enables an indirect observation 
of S continuous in time. The input processes B(t), B^{t) refer to the field before its interaction 
with S, whereas the output processes B out (t) = W (t)B(t)U(t), B out ^(t) = W (t)B^ (t)U(t) 
describe the field after the interaction with S. 

In a balance heterodyne measurement depicted in Fig. 1 the output field, escaping from 
the cavity, is mixed with a strong laser field B\ (t) (local oscillator). We assume that this 
auxiliary field, which does not interact with S, is initially in a coherent state e(f\ ) 2(], 30]. In 
the paper we consider the filtering equation corresponding to the observation of the difference 
of photocurrents generated by the detectors monitoring the fields: 

1 .. .. 1 



B 1 {t) 



V2 



(B° ut (t) + B lo (t)) , B 2 (t) 



V2 



(B out (t)-B lo (t)) 



(9) 



To derive the linear stochastic equation one can use, for example,the Belavkin method of 
generating functional jsj. The generating map, g(k,t), defined by ja |20j 



where 



g(k,t):Z^g(k,t)[Z], 
(V|<7(M)[£]V) = (^®e(f)\G out (k,t)Z^®e(f)) 



t 

G out (M) = (e(/io)|exp{y ek(t')(dA 11 (t')-dA 22 (t': 



|e(/u,)> 



(10) 



IT 



completely determines the observed process up to time t. Here Z t = W(t)ZU(t) is the 
Heisenberg operator of S, ip stands for the initial state of S, k is any integrable c-valued 
function, and e _1 = | /i |- In the limit e — > of a very intense local oscillator field, the 
formula ffTTj) takes the form 



G out (k,t) 



exp{ j k(t')dQ™\t') 



(12) 
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where the output process Q out reads 

Q out (t) = f fe i ^' ) d J B 0Utt (t')+e- i ^' ) d J B 0Ut (t / ) > ), (13) 



<p(t) = arg/io(t). By using (JSJ) and the Ito formula ([6]) one obtains 

dQ out (t) = e^dB^t) + e- [(f,{t) dB{t) + y/jl (e'^aj + e-^ajdt (14) 

and (dQ ont (t)) 2 = dt. An explicit expression for the generating map ( ITUl) one can find by 
solving the differential equation for g(k,t). Using the method described in [^j], one can check 
that the generating map g(k, t) satisfies the equation 



-g(k,t)[Z] = g(k,t)[-(K+^ia^f(t)-^JIaf(t)^Z + ^Za 



-Z{K + v^at/(i) - y/fiaf{t)) + -k\t)Z 

+k{t){yfjia) + J{fy)d*®Z + k(t)Z(^jia + /(t))e- i ^j (15) 

with gj_/c, 0)[Z] = Z and K = j t H + |a^a. The solution to the Eq. ( IT5l) can be written in the 
form [3] 

g(fc,t)[Z] = I 'Gfcq^fflZVfflMQ*), (16) 



where z/ is the probabilistic measure on the set Q consisting of the continuous trajectories 
q = [q(t)\t > 0] of the observed process Q ont (t), restricted to the set f2* = {q l \q G f2} of the 
trajectories q l = [q(r)\r < t] up to t. The stochastic propagator V(t)(q t ) = V(q t ) satisfies 
the stochastic differential equation 

dV(t) = -{K+^Jla ] f{t)+^Jlae- 2i ^ f^V^dt+^Jiae-^V^dQit), 9(0) = I 

(17) 

and 



G(A;,g i ) = G out (k,t)(q) = exp { J k(t')dq(t') 



Hence, the posterior unnormalized wave function ip(t) = V(t)ip of the system S satisfies the 
stochastic dissipative differential equation of the form 



df{t) = -\K + y/Jl (<Jf(t) -af(t)J W(f)d* + V^ae^ (t V(0 d ^(0, ^(°) = V> , 

(19) 

where 

dW(t) = dQ(t) -2Re(e-^ (t) /(t)) dt (20) 
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and W(t) is isomorphic to the standard Wiener process. Let us notice that the integral 
representation (1301) of the generating map, 

(G™\k,t)Z t ) = J Gik.q^Viq^ZViq^Mq 1 ), (21) 

gives for Z = I the mean value of the output process ffl~2]) as the generating function of the 
output probability measure 

dC(?*) = (VtfWVtfWMQ*)- (22) 

Therefore, the posterior wave function ip(t) satisfying the filtering equation ( fl9l) is normal- 
ized to the probability density (V^*)^ |V(g*)^>) of the observed process with respect to the 
probability measure of the input process. From (f2~Tl) we obtain the posterior mean value 
(Z)(q t ) as 

(ZW) = (tptfyZtptf)) , (23) 

where tp(t)(q) = (p(q l ) and <p(t) = (i/j(t)\ip(t))~ 1 ^ 2 i/j(t). For the normalized posterior wave 
function (p(t) we get the nonlinear filtering equation 

d(p{t) = -(^ + ^aV(t)-v^ a M + /^Re 2 ((a) t e- i ^ ) )/2)^(t)dt 
+ / uae-^ (t) Re((a) f e- i ^ ) )^(t)dt+ (^oe"*' - ^/7Re((a) t e" 1 *®)) 



x 



(dW(t) - 2 v ^IRe((a) t e- i ^dt)) , (24) 



where (a) t = (£>(t)|a <p(t)). And if the initial state of S is a mixed one, g{0) = g, then the 
posterior normalized density matrix p[t) satisfies the nonlinear filtering equation of the form 

dp(t) = (~[H,p(t)} - £{a*a,p(t)} + [^Jlaf{t) - y/ji<Sf{t), p(t)} + pap{t)a^j dt 

+ (rfiae-Wfft) + p(t)VJ^^^ m ~ 2 v ^Re((a) t e- i ^)) 

x (dW(t) - 2 v ^IRe((a) t e- i * (i) dt)) , (25) 



where (a) t = Tr[p(t)a]. Eq. (125]) is consistent with the result in [10]. 



4. Posterior evolution of a squeezed coherent state 

We shall show that the coherent state survives the reduction of the state following the 
registered trajectory: the solution to Eq. (JIT?]) for the initial coherent state \a ), a G C, can 
be written in the form \ip{t)) = l(t)\a(t)). Inserting the predicted solution into Eq. ( JT9l and 
making use of the property 

+ , , d\a) 1 . . .„„. 
a j \a) = ^ + 2«l a )' (26) 
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on can write both sides of the equation in terms of linearly independent vectors \a) and -4^-. 
Comparing the coefficients of the corresponding vectors one gets the consistent system of 
differential equations 

da{t) = -(iu+^)a(t)dt- y/jlf(t)dt, (27) 



dl(t) 



l(t) 2 



--dt - ||«(t)| 2 dt - _^a(t)/(t)dt 

a(fjf(t)dt + y/jiaifye-Wt) dW(t), (28) 



for the functions a(t) and l(t) with the initial condition a(0) = «o, /(0) = 1. To solve the 
stochastic equation for l(t) one has to use the Ito rules and the formula 

d, ^-ik dl -m {dl(t))2 - (29) 

Taking into account that (d/(t)) 2 = fxa 2 (t)e~ 21 ^ dt, we obtain the analytical solution of Eq. 
([19]) for the initial coherent state in the form 



\$ao(t)) = exp 



J -f-^l {^W)\ 2 + <t')f{t^dt' 
-V?jT (Re(a(t')TM) + ^a\t')e~^ dt' 

+y/Ji f a(t') e~ im dW(t') \a(t)) , (30) 
Jo J 

where the independent of the noise amplitude a(t) reads 

a(t) = aoe-(^>- y/Jl f e'^)^ f {t')dt' . (31) 



Since any initial state can be represented in the basis of the coherent states 

|^(0)) = - ! d 2 « (ao|^(0))|«o), (32) 



7T _ 

the linearity of Eq. (Tl9|) allows us to write the general solution as 



m)) = - I d 2 a («o|^(0))|^ (t)). 



7T 



The process of measurement changes the state of the system. When the system interacts 
with the external Bose field but the signal is not measured, the state in general becomes 
mixed. The nondemolition observation of the Belavkin type gives the opportunity to get 
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information about the system and allowing to retain some properties of the initial state, for 
example, its purity. 

A coherent state is not the only one state invariant under the filtering equation (fT9l . We 
shall prove that a squeezed coherent state is preserved under the considered observation as 
well. We will discuss the time development of the posterior wave function of S for initial 
state of the form 

|-0(O)> = S(t o )D(a o )\0) = S(£ )\ao) = \£ ,a ) , (34) 

where 

D(oto) = exp [aoo) — o>oa) , «o G C, (35) 

and 

S(£o) = exp Q^a 2 - ^o(a f ) 2 ) , £o = ^° Qo e C. (36) 

The state is generated by displacing the vacuum state and then by squeezing. The amount of 
squeeze is described by p = |£ | which is called the squeeze factor. Some details of description 
and detection of squeezed states of light one can find, for instance, in 11). 

The method of computing a posterior dynamics for the initial squeezed coherent state by 
making a simple ansatz like in the previous case is cumbersome and laborious. Therefore, we 
shall use a more efficient method which allows to avoid arduous computation of a stochastic 
phase of posterior state. 

We make use of the eigenvalue equation 

S(t)atf(£)\£,a) = a\£,a), (37) 

which can be readily. fond from the definition (134)) . The operator expansion theorem allows 
one to check that 11] 

S(£)a&(£) = aT 1 + a^T 2 , (38) 

where Ti = cosh g, T 2 = e ie sinh g. Let us notice that if the system S remains in the squeezed 
coherent state at any time instant t > 0, then 

S(Z(t))a£fi(Z(t))$(t) = [ar^ + atr^t)]^, (39) 

and 

S(£(t + dt))aS\£(t + dt))$(t + dt) 

= [aT 1 (t + dt) + a j T 2 (t + dt)] ${t + dt), (40) 

where \ip(t)) = l(t)S(£(t))\a(t)), have to be fulfilled. Eqs. (1391 and (1401 can be reduced to a 
single condition of the form 

[a(r 1 (t) + dr 1 (t)) + a t (T 2 (t)+dT 2 (t))-a(t)-da(t)] d$(t) 
+ (adT 1 (t) + a^dT 2 (t)-da(t))$(t) = 0. (41) 
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Finally inserting of the increment dip(t) given by Eq. (fl9|) into Eq. (I4ip we obtain the set of 
the two differential equations: 



r 2 (0 [-r\(0 (iw + 1) dt +/ie- 2i ^r 2 (t)dt +dr x (t) 



o 



-ri(t) [r 2 (t) + dt + dr 2 (t) 



(42) 



a(t)ri(t) [-rx(t) dt+//e- 2i ^r 2 (t)dt+dr!(t)" 
-a(t)Ttfj[r 2 (t) dt+dr 2 (t)] -^r!(t)/(t)di 

- v ^r 2 (t)J(t)dt - v//ir 2 (t)e- i<Wt) dW(t) - da(i) = 



(43) 



with the initial condition: 1^(0) = cosh £ , T 2 (0) = e l9 ° sinh g , a(0) = a . The last step 
requires left-multiplying the Eq. (JS) by S f (£(t)) and use of the transformation 



tf{£)aS(£) = aT 1 -a^T 2 



(44) 



which one can easily get from (1381) . The Eqs. (|42p . (j43j) form a consistent set of equations 
and this is completes the proof. 

Eq. ( l42|) can be rewritten in terms of the function r(£) = T 2 (t)/ri(t) as 



A r( t) = -2(ic + |)r(t) + / ,e- 2i ^r 2 (t). 

The general solution to the Riccati differential equation 045)) reads 

r(o)e~( 2itj+ ^ 4 



r(t) 



l-/ir(0)/e^( 2i ^)''- 2i ^')dt' 



(45) 



(46) 



In particular, for the phase <p(t) = ir/2 — oj t, one gets 

{2iu - 2iw + A*)r(0) 



T(t) 



e(2i^+M)t[2ia; _ 2iw + ^(1 + r(0))] - ,ur(0)e 2i "°* 
The integration of (143 p yields the stochastic amplitude 

1 

a(t) = 



(47) 



a \ 1 - T 



x exp ^— io^t -^t + fij* e- 2i * (i,) r(t')dt^ 

- ^ jT exp f- (iu; + |) (t - s) + fi jf e- 2i ^T(t')dt' 
x (f(s)ds + T{s)J{s)ds + e-^ {s) T{s)dW{s 



(48) 



Therefore the posterior mean values of the optical quadratures X = (a + a')/2 and Y 
(a — a ! )/2i for the initial squeezed coherent state given as 



(Y)t 



Re la(t) - a(t)T(t) 



Im la(t) -a(t)T(t) 



(49) 
(50) 

(51) 
(52) 

(53) 



x/i-irWI 2 

depend on the measurement noise, whereas the uncertainties of quadratures 

AX(t) = (4:ReK(t)y 1/2 , 
AY{t) = |«(t)| (4Re/t(t))" 1/2 , 

where 

«f« " 1 + r(t) 
K{t) ~ l-Y{t) 

remain deterministic. Moreover, the expressions for the mean value of optical quadratures 
include the parameter of the initial coherent state of Bose field, whereas the formulae for 
the uncertainties are exactly the same as for the case when the Bose field is initially in the 
vacuum state 121 ] . 

The time dependence of the uncertainties AX and AY has been illustrated by the para- 
metric plots presented in Fig. 2. They show the dynamics of AX and AY as functions of 
the dimensionless time r = out (0 < r < 100), for ujq = 0.05, T(0) = 0.8 and three values 
of [i. The uncertainties oscillate approaching the asymptotic values (AX = AY = |). The 
lager parameter \x the asymptotic values are reached faster. 

In order to study the probability density of the output process (1221 one has to find the 
norm ||^(^)|| = I KOI- O ne can check that 



dl*(*)l 5 

W)\ 2 









a(t)r(t)) 


e -i<f>(t) 




r(*)| 2 



dW(t) 



(54) 



Hence 



\Kt)\'' 



rt Re 

exp ^ lyfJL J 



a{t') -a(t')T{t') e"^*') 



dW(t') 



(55) 



-2fi 



t Re 2 



a(t') - a(t')r(t') e- 1 ^*') 



|r(f)| 5 



dt' 



Let us stress that the formula (1221) gives the probability of observed result q 1 £ fi*. 
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We should remark that, due to lim^oo T(t) = 0, the system asymptotically approaches 
the coherent state with the amplitude 



Hence, any memory about the initial condition is lost after a transitional period. By (|33|) 
the asymptotic posterior state for any initial state is the coherent state with the amplitude 



5. Final remarks 

We have presented the derivation of the Belavkin filtering equation for a single-mode field 
in a cavity interacting with the Bose field (measuring apparatus) initially prepared in a 
coherent state. The considered balanced heterodyne observation is intense and is therefore 
considered as a diffusion one. In contrast to the case of the measuring apparatus prepared 
initially in the vacuum state, when the initial squeezed coherent state is preserved [12[] and 
driven to the vacuum asymptotic state, in the case studied in this paper the initial squeezed 
coherent state is asymptotically driven to the coherent one. This asymptotic state does not 
depend on the initial parameters of the initial state of the field in the cavity. Moreover, we 
have proved that any initial state described by a square integrable wave function relaxes to 
the coherent state, with the amplitude dependent on the coupling constant and the initial 
coherent state of the apparatus. The driving force from the field can control the system and 
drive its state to the coherent one with the given amplitude. 

References 

[1] R. L. Hudson and K. R. Parthasarathy, "Quantum Ito's formula and stochastic evolutions, 



Commun. Math. Phys. 93, 301-323 (1984). 

[2] K. R. Parthasarathy, An Introduction to Quantum Stochastic Calculus (Birkhauser, 1992). 

[3] V. P. Belavkin, "A posterior Schrodinger equation for continuous nondemolition measure- 
ment", J. Math. Phys. 31, 2930-2934 (1990). 

[4] A. Barchielli and V. Belavkin, "Measurement continuous in time and a posteriori states 
in quantum mechanics", J. Phys. A: Math. Gen. 24, 1495-1514 (1991). 

[5] V. P. Belavkin, "Measurement, filtering and control in quantum open dynamical sys- 
tems", Rep. Math. Phys. 43, 405-425 (1999). 

[6] V. P. Belavkin, "Quantum causality, decoherence, trajectories and information" , Rep. 
Prog. Phys. 65, 353-420 (2002). 




(56) 



(EH}. 



li 



[7] A. Barchielli, "Direct and heterodyne detection and other applications of quantum 
stochastic calculus to quantum optics", Quantum Opt. 2, 423-441 (1990). 

[8] A. Barchielli and A. M. Paganoni, "Detection theory in quantum optics: stochastic rep- 
resentation", Quantum Opt. 8, 133-156 (1996). 

[9] J. E. Gough and C. Kostler, "Quantum filtering in coherent states", Commun. Stoch. 
Anal. 4, 505-521 (2010). 



[10 

[11 
[12 

[13 

[14 



[15 

[16 
[17 

[18 
[19 
[20 



A. Dabrowska and P. Staszewski, "Filtering equation for measurement of a coherent 
channel", J. Opt. Soc. Am. B 28, 1238-1244 (2011). 

Ch. Gerry and P. Knight, Introductory Quantum Optics (Cambridge University Press, 
2005). 

A. Dabrowska and P. Staszewski, "Squeezed coherent state undergoing a continuous 
nondemolition observation", Phys. Lett. A 375, 3950-3955 (2011). 

H. Carmichael, An Open Systems Approach to Quantum Optics ( Springer- Verlag Berlin 
Heidelberg, 1993). 

A. C. Doherty and H. Mabuchi, "Atoms in micro cavities: quantum electrodynamics, 
quantum statistical mechanics, and quantum information science" , in Optical microcavi- 
ties, K. Vahala, ed. (World Scientific Press, Singapore, 2004), pp. 367-414. 

H. Wiseman and G. J. Milburn, Quantum measurement and control (Cambridge Uni- 
versity Press, 2010). 

T. Brun, "A single model of quantum trajectories", Am. J. Phys. 70, 719-737 (2002). 

R. van Handel, J. K. Stockton, and H. Mabuchi, "Modelling and feedback control 
design for quantum state preparation", J. Opt. B: Quantum Semiclass. Opt. 7, S179- 
S197 (2005). 

K. Jacobs and D. Steck, "A straightforward introduction to continuous quantum 
measurement", Contemporary Physics 47, 279-303 (2006). 

L. Bouten, R. van Handel, and M. James, "An introduction to quantum filtering", SI AM 
J. Control Optim. 46, 2199-2241 (2007). 

A. Barchielli, "Continual measurements in quantum mechanics and quantum stochastic 
calculus", in Open Quantum Systems III. Recent Developments, S. Attal, A. Joye, and 
C.-A. Pillet, eds. (Berlin Heidelberg: Springer LNM, 2006), pp. 207-288. 



12 



[21] V. P. Belavkin and M. Guta, eds. Quantum stochastics and information, (World Scien- 
tific, 2006). 

[22] A. Barchielli and M. Gregoratti, Quantum Trajectories and Measurement in Continuous 
Time: The Diffusive Case, ( Springer- Verlag Berlin Heidelberg, 2009). 

[23] C. W. Gardiner and P. Zoller, Quantum noise (Springer- Verlag Berlin Heidelberg, 2010). 

[24] H. Mabuchi, J. Ye, and H. J. Kimble, "Full observation of single-atom dynamics in 
cavity QED", Appl. Phys. B 68, 1095-1108 (1999). 

[25] M. A. Armen, J. K. Au, J. K. Stockton, A. C. Doherty, and H. Mabuchi, 'Adaptive 
homodyne measurement of optical phase", Phys. Rev. Lett. 89, 133602-133602 (2002). 

[26] J. M. Geremia, J. K. Stockton, A. C. Doherty, and H. Mabuchi, "Quantum Kalman 
Filtering and the Heisenberg Limit in Atomic Magnetometry" , Phys. Rev. Lett. 91, 
250801-250804 (2003). 

[27] J. M. Geremia, J. K. Stockton, and H. Mabuchi, "Tensor polarizability and dispersive 
quantum measurement of multilevel atoms", Phys. Rev. A 73, 042112-042126 (2006). 

[28] P. Goetsch and R. Graham, "Linear stochastic wave equation for continuously measured 
quantum system", Phys. Rev. A 50, 5242-5255 (1994). 

[29] M. J. Collet and C. W. Gardiner, "Input and output in damped quantum systems: Quan- 
tum stochastic differential equations and the master equation", Phys. Rev. A 31, 3761- 
3774 (1985). 

[30] A. Barchielli, "Direct and heterodyne detection and other applications of quantum 
stochastic calculus to quantum optics", Quantum Opt. 2, 423-441 (1990). 



13 



photodetector 2 



system S 



output field 



beam 
splitter 



photodetector 1 



local 
oscillator 



Fig. 1. Balanced heterodyne detection 
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Fig. 2. The dependence of AX and AY on the dimensionless time r = out is 
displayed for ojq = 0.05, T(0) = 0.8 and for three values of fi: 0.01 (a), 0.04 
(b), and 0.08 (c). 
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